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Abstract — Lorentz Invariance (LI) is the founding postulate 
of Einstein's 1905 theory of relativity, and therefore at the 
heart of all accepted theories of physics. It characterizes the 
invariance of the laws of physics in inertial frames under changes 
of velocity or orientation. This central role, and indications 
from unification theories [1], [2], [3] hinting toward a possible 
LI violation, have motivated tremendous experimental efforts 
to test LI. A comprehensive theoretical framework to describe 
violations of LI has been developed over the last decade [4]: the 
Lorentz violating Standard Model Extension (SME). It allows 
a characterization of LI violations in all fields of present day 
physics using a large (but finite) set of parameters which are all 
zero when LI is satisfied. All classical tests (e.g. Michelson-Morley 
or Kennedy-Thorndike experiments [5], [6]) can be analyzed 
in the SME, but it also allows the conception of new types 
of experiments, not thought of previously. We have carried out 
such a conceptually new LI test, by comparing particular atomic 
transitions (particular orientations of the involved nuclear spins) 
in the 133 Cs atom using a cold atomic fountain clock. This allows 
us to test LI in a previously largely unexplored region of the 
SME parameter space, corresponding to first measurements of 
four proton parameters and an improvement by 11 and 12 orders 
of magnitude on the determination of four others. In spite of the 
attained accuracies, and of having extended the search into a new 
region of the SME, we still find no indication of LI violation. 

The experiment reported here tests for a violation of LI 
by searching for a variation of atomic transition frequen- 
cies as a function of the orientation of the spin of the 
involved atomic states (clock comparison or Hughes-Drever 
experiment). Based on the SME analysis of numerous atomic 
transitions in [7], [8] we have chosen the measurement of a 
particular combination of transitions in the 133 Cs atom. This 
provides good sensitivity to the quadrupolar SME energy shift 
of the proton (as defined in [7]) in the spin | Cs nucleus, 
whilst being largely insensitive to magnetic perturbations (first 
order Zeeman effect). The corresponding region of the SME 
parameter space has been largely unexplored previously, with 
first limits on some parameters set only very recently [9] by 
a re-analysis of the Doppler spectroscopy experiment (Ives- 
Stilwell experiment) of [10]. Given the large improvement 
(11 and 12 orders of magnitude on four parameters) we 
obtain with respect to those results, and the qualitatively 
new region explored (first measurements of four parameters), 
the experiment had comparatively high probability for the 
detection of a positive Lorentz violating signal. However, our 
results clearly exclude that possibility. 



The SME was developed relatively recently by Kostelecky 
and co-workers [4], motivated initially by possible Lorentz 
violating phenomenological effects of string theory [1]. It 
consists of a parameterized version of the standard model 
Lagrangian that includes all Lorentz violating terms that can 
be formed from known fields, and includes (in its most recent 
version) gravity. In its minimal form the SME characterizes 
Lorentz violation in the matter sector by 44 parameters per 
particle [7], [8], of which 40 are accessible to terrestrial 
experiments at first order in u©/c [8] (u® is the orbital 
velocity of the Earth and c = 299792458 m/s). All SME 
parameters vanish when LI is verified. Existing bounds for the 
proton (p + ), neutron (n) and electron (e~) come from clock 
comparison and magnetometer experiments using different 
atomic species (see [7] and references therein, [11], [12], [13], 
[14]), from resonator experiments [15], [6], [16], and from 
Ives-Stilwell experiments [9], [10]. They are summarized in 
tab. U together with the results reported in this work. 

For our experiment we use one of the laser cooled fountain 
clocks operated at the Paris observatory, the 133 Cs and 87 Rb 
double fountain F02 (see [17] for a detailed description). We 
run it in Cs mode on the \F = 3} <-> \F = 4) hyperfine 
transition of the 66*1/2 ground state. Both hyperfine states are 
split into Zeeman substates rap = [—3,3] and rap = [—4,4] 
respectively. The clock transition used in routine operation is 
\F = 3,m F = 0) \F = 4,m F = 0) at 9.2 GHz, which 
is magnetic field independent to first order. The first order 
magnetic field dependent Zeeman transitions (|3, i) <-» |4, z) 
with i = ±1, ±2, ±3) are used regularly for measurement and 
characterization of the magnetic field, necessary to correct the 
second order Zeeman effect of the clock transition. In routine 
operation the clock transition frequency stability of F02 is 
1.6 x lO-^r" 1 / 2 , and its accuracy 7 x lO" 16 [17], [18]. 

A detailed description of the minimal SME as applied to the 
perturbation of atomic energy levels and transition frequencies 
can be found in [7], [8]. Based on the Schmidt nuclear model 1 
one can derive the SME frequency shift of a Cs \3,m F ) <-> 

1 As discussed in [7] the Schmidt nuclear model only allows an approximate 
calculation of the SME frequency shift, with more complex models generally 
leading to dependences on additional SME parameters. Nonetheless, the model 
is sufficient to derive the leading order terms and has been used for the analysis 
of most experiments providing the bounds in Tab. U 



TABLE I 

Orders of magnitude of present limits (in GeV) on Lorentz violating parameters in the minimal SME matter sector and 
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for the quantization magnetic field in the negative z direction 
(vertically downward) in the lab frame. The first three terms 
in Q are Lorentz violating SME frequency shifts, the last two 
describe the first and second order Zeeman frequency shift (we 
neglect B 3 and higher order terms). The tilde quantities are 
linear combinations of the SME matter sector parameters of 
tab. U in the lab frame, with p,e standing for the proton and 
electron respectively. The quantities f3 w ,5 w ,K W) / y w , X w (with 
w = p, e) depend on the nuclear and electronic structure, they 
are given in tab. II of [8], h is Planck's constant, B is the 
magnetic field seen by the atom, — 7.008 Hz nT" 1 is the 
first order Zeeman coefficient, = 427.45 x 10 s Hz T~ 2 
is the second order coefficient [19]. The tilde quantities in 
Q are time varying due to the motion of the lab frame 
(and hence the quantization field) in a cosmological frame, 
inducing modulations of the frequency shift at sidereal and 
semi-sidereal frequencies, which can be searched for. 

From we note that the nip — clock transition is 
insensitive to Lorentz violation or the first order Zeeman 
shift, while the Zeeman transitions (mp ^ 0) are sensitive 
to both. Hence, a direct measurement of a Zeeman transition 
with respect to the clock transition allows a LI test. However, 
such an experiment would be severely limited by the strong 
dependence of the Zeeman transition frequency on B, and its 
diurnal and semi-diurnal variations. To avoid such a limitation, 
we "simultaneously" (see below) measure the mp = 3, mp = 
—3 and mp = transitions and form the combined observable 
v c = 2/ + 3 + i^_3 — 2^o. From Q 



v c = ^K p cP - -KfB* (2) 

where we have used j p — —K p /9 from [8] (K v ~ 10~ 2 in the 
Schmidt nuclear model). This observable is insensitive to the 
first order Zeeman shift, and should be zero up to the second 
order Zeeman correction and a possible Lorentz violating shift 
in the first term of 0. 

The lab frame parameter cP q can be related to the con- 
ventional sun-centered frame parameters of the SME (the 
parameters of tab. |l} by a time dependent boost and rotation 
(see [8] for details). This leads to a general expression for the 
observable v c of the form 

v c = A + C^ e cos(w e T e ) + 5 w@ sin(w e T e ) (3) 

+ C 2we cos(2a> e T e ) + S 2we sin(2w e T e ), 

where is the frequency of rotation of the Earth, T ffi is time 
since 30 March 2005 1 lh 19min 25s UTC (consistent with the 
definitions in [20]), and with A, C W{B , S u@ , C 2uJm and S 2u , m 
given in tab.||T|as functions of the sun frame SME parameters. 
A least squares fit of Q to our data provides the measured 
values given in tab. |ll] and the corresponding determination of 
the SME parameters. 

The F02 setup is sketched in fig-d Cs atoms effusing from 
an oven are slowed using a chirped counter propagating laser 
beam and captured in a lin _L lin optical molasses. Atoms are 
cooled by six laser beams supplied by pre adjusted optical fiber 
couplers precisely attached to the vacuum tank and aligned 
along the axes of a 3 dimensional coordinate system, where 
the (111) direction is vertical. Compared to typical F02 clock 
operation [17], the number of atoms loaded in the optical 
molasses has been reduced to 2 x 10 7 atoms captured in 30 ms. 
This reduces the collisional frequency shift of v c to below 
0.1 mHz, and even less for its variation at w® and 2W®. 

Atoms are launched upwards at 3.94 m.s -1 by using a 
moving optical molasses and cooled to ~ 1 //K in the moving 



TABLE II 

Coefficients of (5J to first order in p = v 9 /c, where n$ is the angular frequency of the Earth' s orbital motion, T is time since 
the March equinox, \ = 41.2° is the colatitude of our lab, and rj = 23.3° is the inclination of the Earth's orbit. The measured 

VALUES (IN mHz) ARE SHOWN TOGETHER WITH THE STATISTICAL (FIRST BRACKET) AND SYSTEMATIC (SECOND BRACKET) UNCERTAINTIES . FOR OUR 
RELATIVELY SHORT DATA SET (21 DAYS) WE NEGLECT THE SLOW VARIATION DUE TO THE ANNUAL TERMS AND TAKE ~ 0.34 rad. 
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frame by adiabatically decreasing the laser intensity and in- 
creasing the laser detuning. Atoms are then selected by means 
of a microwave excitation in the selection cavity performed 
in a bias magnetic field of ~ 20 /iT, and of a push laser 
beam. Any of the 1 3, Trip) states can be prepared with a high 
degree of purity (few 1CP 3 ) by tuning the selection microwave 
frequency. 52 cm above the capture zone, a cylindrical copper 
cavity (TEon mode) is used to probe the |3, mp) <-> |4, mj?) 
hyperfine transition at 9.2 GHz. The Ramsey interrogation 
method is performed by letting the atomic cloud interact with 
the microwave field a first time on the way up and a second 
time on the way down. After the interrogation, the populations 
Np = 4 and Np = % of the two hyperfine levels are measured 
by laser induced fluorescence, leading to a determination of 
the transition probability P = ^/(N^ + N4) which is in- 
sensitive to atom number fluctuations. One complete fountain 
cycle from capture to detection lasts 1045 ms in the present 
experiment. From the transition probability, measured on both 
sides of the central Ramsey fringe, we compute an error signal 
to lock the microwave interrogation frequency to the atomic 
transition using a digital servo loop. The frequency corrections 
are applied to a computer controlled high resolution DDS 
synthesizer in the microwave generator. These corrections are 
used to measure the atomic transition frequency with respect 
to the local reference signal used to synthesize the microwave 
frequency. 

The homogeneity and the stability of the magnetic field in 
the interrogation region is a crucial point for the experiment. 
A magnetic field of 203 nT is produced by a main solenoid 
(length 815 mm, diameter 220 mm) and a set of 4 compen- 
sation coils. These coils are surrounded by a first layer of 3 
cylindrical magnetic shields. A second layer is composed of 2 
magnetic shields surrounding the entire experiment (optical 
molasses and detection zone included). Between the two 
layers, the magnetic field fluctuations are sensed with a flux- 
gate magnetometer and stabilized by acting on 4 hexagonal 
coils. The magnetic field in the interrogation region is probed 
using the |3, 1) <-> |4, 1) atomic transition. Measurements of 
the transition frequency as a function of the launch height show 
a peak to peak spatial dependence of 230 pT over a range of 
320 mm above the interrogation cavity with a variation of 
< 0.1 pT/mm around the apogee of the atomic trajectories. 



Measurements of the same transition as a function of time at 
the launch height of 791 mm show a magnetic field instability 
near 2 pT for an integration time of r =1 s. The long term 
behavior exhibits residual variations of the magnetic field (~ 
0.7 pT at t =10000 s) induced by temperature fluctuations 
which could cause variations of the current flowing through 
solenoid, of the solenoid geometry, of residual thermoelectric 
currents, of the magnetic shield permeability, etc... 

The experimental sequence is tailored to circumvent the 
limitation that the long term magnetic field fluctuations could 
cause. First |3, — 3) atoms are selected and the 3,— 3} <-> 

4, —3) transition is probed at half maximum on the red 
side of the resonance (0.528 Hz below the resonance center). 
The next fountain cycle, |3, +3) atoms are selected and the 

3, +3) <-> 4, +3) transition is also probed at half maximum 
on the red side of the resonance. The third and fourth fountain 
cycles, the same two transitions are probed on the blue side of 
the resonances (0.528 Hz above the resonance centers). This 
4180 ms long sequence is repeated so as to implement two 
interleaved digital servo loops finding the line centers of both 
the |3, —3) <-> |4, -3} and the |3,+3) «-> |4, +3) transitions. 
Every 400 fountain cycles, the above sequence is interrupted 
and the regular clock transition |3, 0} <-> |4, 0) is measured for 
10 s allowing for an absolute calibration of the local frequency 
reference with a suitable statistical uncertainty. Using this 
sequence, magnetic field fluctuations over timescales > 4 s 
are rejected in the combined observable v c and the stability 
is dominated by the short term (r < 4 s) magnetic field 
fluctuations. 

We have taken data implementing the experimental se- 
quence described above over a period of 21 days starting on 
march 30, 2005. The complete raw data (no post-treatment) is 
shown in fig.|2] each point representing a ^432 s measurement 
sequence of + ^-3 — 2^o as described above. The inset 
in fig. [2] shows the frequency stability of the last continuous 
stretch of data (^10 days). We note the essentially white 
noise behavior of the data on that figure, indicating that 
the experimental sequence successfully rejects all long term 
variations of the magnetic field or of other perturbing effects. 
A least squares fit of the model to the complete data 
provides the 5 coefficients and statistical uncertainties given 
in tab. HD 



We note a statistically significant offset of the data from 
zero (—5.7(0.05) mHz). This is partly due to the second 
order Zeeman shift (second term in (|2ji) which amounts to 
—2.0 mHz. The remaining —3.7 mHz are either due to 
a systematic effect or indicate a genuine Lorentz violating 
signal. 

The dominant systematic effect in our experiment is most 
likely a residual first order Zeeman shift. This arises when 
the nip = —3 and mp = +3 atoms have slightly different 
trajectories in the presence of a magnetic field gradient. The 
result is a difference in first order Zeeman shift and hence 
incomplete cancelation in the combined observable v c . To 
test this hypothesis we have taken 3 hours of time of flight 
(TOF) measurements of the atoms in the different Zeeman 
states. We observe a significant TOF difference of 155 /xs 
(when launching at 4.25 m/s) between the nip = —3 and 
nip = +3 states, most likely due to slight magnetic and optical 
differences in the laser trapping and cooling during the final 
adiabatic phase. 

We model a residual first order Zeeman effect using a Monte 
Carlo (MC) simulation with the measured vertical magnetic 
field map, and assuming that the complete TOF difference is 
due to an initial vertical spatial offset (~ 660 /im) between the 
atoms rather than an offset in launch velocity. In this scenario, 
which maximizes the resulting frequency shift, we obtain a 
frequency offset of only ~ 5 [jHz in v c . 

Another indication for magnetic field inhomogeneities 
comes from measurements of the contrast of the Ramsey 
fringes for the different Zeeman transitions. The measured 
contrasts are 0.94, 0.93, 0.87, and 0.75 for the m f = 
0, mp — +1, mp = +2, and mp = +3 transitions 
respectively. The observed differences cannot be explained 
by the known vertical magnetic field gradient alone, but 
are due to an additional horizontal field gradient, probably 
caused by residual magnetization inside the magnetic shields 
or lack of symmetry in the shielding or compensation. Our 
MC simulation reproduces the measured contrasts when we 
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include a constant horizontal gradient of ~ 6 pT/mm. To 
check this value we have slightly tilted the fountain (by up to 
~ 1800 /irad off the vertical) and measured the corresponding 
change in frequency of the mp = +1 Zeeman transition. The 
obtained field gradient is somewhat smaller (~ 2 pT/mm) but 
we consider the latter method less accurate, and conservatively 
use the larger value obtained from the contrast measurements. 

Including the thus determined horizontal field gradient in 
a complete MC simulation we obtain a total offset in v c of 
~ 26 mHz, assuming that the horizontal spatial offset between 
the mp = —3 and mp = +3 atoms is identical to the 
measured vertical offset (660 /im). We consider this to be an 
upper limit (the horizontal separation is likely to be less than 
the vertical one due to the absence of gravity), and take it as 
the systematic uncertainty on the determined offset (A in tab. 

m. 

To determine the systematic uncertainty on the sidereal 
and semi-sidereal modulations of v c we have measured the 
variation of the mp = TOF at those frequencies and 
taken the result as the maximum variation of the mp = +3 
vs mp = —3 TOF difference (note that the real value is 
probably smaller, due to common mode variations). The best 
fit amplitudes of the sidereal and semi-sidereal TOF variation 
are 3.0 /is and 2.9 /.is respectively, which leads to an upper 
limit of 0.35 mHz on the systematic uncertainties of C u@ , 



tab. ED 



Finally we use the five measurements and the relations in 
tab. [n] to determine the values of the eight SME parameters. 
In doing so, we assume that there is no correlation between 
the three cpj parameters and the other five parameters, and 
determine them independently. The results are given in tab. 
EH 

In conclusion, we have carried out a test of Lorentz invari- 
ance in the matter sector of the minimal SME using Zeeman 
transition in a cold 133 Cs atomic fountain clock. We see 
no indication of a violation of LI at the present level of 
experimental uncertainty. Using a particular combination of 
the different atomic transitions we have set first limits on four 
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Fig. 1. Schematic view of an atomic fountain. 



00 20x10* 4,0x10* 6.0x10* 8.0x10* 1.0x10* 1.2x10* 1.4x10* 1.6x10* 1.8x10* 
time (s) 



Fig. 2. Raw data of the measurements of (f+3 + v~z — 2fo) spanning ~ 21 
days. The inset shows the stability (Alan deviation) of the last continuous 
stretch of data (~ 10 days). 



TABLE III 

Results for SME Lorentz violating parameters c for the 

PROTON, IN GeV AND WITH J = X,Y, Z. 



c Q = -0.3(2.2) xl(T^ 
cj = (0.06(0.7), -0.2(0.7), 
c TJ = (0.1(1.3), -0.04(1.7). 



c_ = -0.2(1.6) xlO" 24 
0.1(1.6)) xlO" 24 
-0.1(1.0)) xlQ- 20 



proton SME parameters and improved previous limits [9] on 
four others by 11 and 12 orders of magnitude. 

Continuing the experiment regularly over a year or more 
will allow statistical decorrelation of the three ctj parameters 
from the other five, due to their modulation at the annual 
frequency (fi©T terms in tab. HJ. Further improvements, and 
new measurements, could come from the unique capability of 
our fountain clock to run on 87 Rb and 133 Cs simultaneously. 
The different sensitivity of the two atomic species to Lorentz 
violation (see [8]) and magnetic fields, should allow a mea- 
surement of all SME parameters in Q in spite of the presence 
of the first order Zeeman effect. Ultimately, space clocks, like 
the planned ACES mission [21] will provide the possibility 
of carrying out similar experiments but with faster (90 min 
orbital period) modulation of the putative Lorentz violating 
signal, and correspondingly faster data integration. 
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